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Abstract 

We have studied dynamical properties and quantum tunneling in asymmetric double- well (DW) 
systems, solving Schrodinger equation with the use of two kinds of spectral methods for a given 
initial Gaussian wavepacket. Time dependences of wavefunction, averages of position and momen- 
tum, the auto-correlation function, an uncertainty product and the tunneling probability have been 
calculated. Our calculations have shown that (i) the tunneling probability is considerably reduced 
by a potential asymmetry AJ7, (ii) a resonant tunneling at | AJ7| ~ k huo is not realized for a mo- 
tion starting from the lower minimum of potential wells (k = 0, 1, 2, • • • ; uj\ oscillator frequency at 
minima), (iii) the reduction of the tunneling probability by an asymmetry is less significant for an 
initial Gaussian wavepacket with narrower width, and (iv) the uncertainty product (5x 2 )(5p 2 ) in 
the resonant tunneling state is larger than that in the non-resonant tunneling state. The item (i) 
is consistent with previous studies and the items (ii) — (iv) are new findings in the present study. 
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I. INTRODUCTION 



Double-well (DW) systems have been extensively studied in a wide range of fields in- 
cluding physics, chemistry and biology (for a recent review on DW systems, see Ref. [1]). 
Quantum tunneling is one of the most fascinating phenomenon in DW systems [2]. Much 
experimental and theoretical studies have been made in tunneling of a quantum particle in 
DW systems. Quantum tunneling of a particle is possible from one-side well to the other- 
side well through classically forbidden region. Well-known old examples of DW systems 
include an inversion of anmonia molecule. In recent years, there has been an advance in 
the experimental study on macroscopic quantum tunneling such as Josephson junction and 
Bose-Einstein condensation in a double trap. 

DW potential needs not be symmetric and it may be asymmetric in general. In many 
experiments, the asymmetry of the DW potential can be changed by modifying external 
parameters. However, most theoretical studies have been made for symmetric DW systems, 
and asymmetric systems have received less theoretical attention than symmetric ones [SHE] • 
This is because solving an asymmetric DW system is more difficult than a symmetric one. 
Theoretical studies on asymmetric DW systems have been made based on various approx- 
imate methods like the WKB for simplified artificial DW potentials which are analytically 
tractable but not realistic [2]. By using such DW potentials, Weiner and Tse [3], and Nieto 
et al [4J showed that although the tunneling probability is significantly reduced by the 
potential asymmetry, it is enhanced when the asymmetry meets the resonance condition. 
Mugai et al [5] studied the fractal nature of the trajectory in asymmetric DW systems. 
By using WKB, Song [6j studied an asymmetric DW system where the difference of the 
potential minima is close of a multiple of huj (harmonic frequency in the wells). Rastelli [7] 
obtained a semi-classical formula for the tunneling amplitude in asymmetric DW systems 
with the use of WKB method. Conventional theories for DW systems have adopted the 
two-level approximation where the initial state in one-dimensional system is assumed to be 
given by ^(x, 0) = [^o(^) ~~ ^l^jj/v^j denoting the vth {y = 0, 1) eigenfunction. In 

order to discuss the tunneling probability in asymmetric DW systems, Cordes and Das [8] 
proposed a generalized two-level approximation: related discussion will be given in Sec. III. 

For a study on dynamics of wavepacket or tunneling in DW systems, it is necessary to 
solve the time-dependent Schodinger equation subject to appropriate initial and boundary 
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conditions [9]. In the past when quantum mechanics was born, it was very difficult to numer- 
ically solve the time-dependent Schodinger equation even for a simple potential except for a 
harmonic oscillator (HO) potential. We had to develop approximation methods applicable 
to simple tractable DW models although they are not necessarily realistic. In recent years, 
however, there has been significant development in computer and its software. It is now 
possible for us to solve the time-dependent Schrodinger equation with a sufficient accuracy, 
by using convenient packages such as MATHEMATICA, MATLAB and Maple. 

The purpose of the present study is to numerically study dynamics of Gaussian wavepack- 
ets and to examine the effect of the asymmetry on quantum tunneling in asymmetric DW 
systems. Quite recently it has been pointed out that a potential asymmetry of a DW system 
has significant effects on its specific heat [10]. We expect that it is the case also on dynam- 
ical properties of DW systems. We will solve the time-dependent Schodinger equation by 
the spectral method for a given initial squeezed Gaussian wavepacket [HJ [12] , adopting the 
realistic quartic DW potential. In order to investigate the influence of the initial state on 
dynamical properties, we adopt two squeezed Gaussian wavepackets with different parame- 
ters. 

The paper is organized as follows. In Sec. II, we will briefly mention the calculation 
method employed in our study [13]. In solving the time-dependent Schrodinger equation, 
we have adopted the two kinds of spectral method A [Eq. d7|)] and spectral method B 



[Eq. (13)] with energy matrix elements evaluated for a finite size N m (= 30). By using the 
spectral method A, we have calculated time-dependences of the magnitude of wavefunction, 
expectation values of position and momentum, the auto-correlation function, the uncertainty 
product and the tunneling probability, whose results are reported in Sec. III. In Sec. IV the 
tunneling probability is discussed with the use of the spectral method B. We discuss also 
wavepacket dynamics when an initial Gaussian wavepacket starts from the top of the DW 
potential. Sec. V is devoted to our conclusion. 
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II. THE ADOPTED CALCULATION METHOD 



A. Spectral method A 

We consider a quantum system whose Hamiltonian is given by 

+ (1) 

where m, x and p express a mass, position and momentum, respectively, of a particle, 
and U(x) denotes the DW potential. Various approximate analytical and numerical meth- 
ods have been proposed to solve the Schrodinger equation [9j. In our study we solve the 
Schrodinger equation, treating the Hamiltonian H as 

H = H + U(x) - U (x) = H + V(x), (2) 

with 

2 

#o = ^- + U (x), (3) 
zrn 

V{x) = U{x) - U (x). (4) 

Here Hq denotes the Hamiltonian with an appropriate potential Uq : for which eigenfunction 
and eigenvalue are assumed to be given by (j) n (x) and E$ ni respectively. 

For the stationary state, we expand the eigenfunction iff(x) in terms of (j) n {x) 

V(x) = ^c n n (x), (5) 

n=0 

leading to the secular equation 

Ec m ^ ^ H rnn c n) (6) 

n=0 

where E denotes the eigenvalue and N m is the maximum quantum number. 

For the time-dependent state, we adopt the spectral method A, in which the eigenfunction 
iff(x,t) is expanded in terms of (j) n {x) 

= ^c n (t)0 n (x), (7) 

n=0 
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where c n {t) stands for the time-dependent expansion coefficient obeying equations of motion 
given by 

ih ojn^ = Hmn Cn{ t) ( m = to N m ), (8) 

n=0 

with 

/oo 
0m(^)* V(x) <j) n (x) dx. (9) 
-co 

Equation ^ expresses the (N m + 1) first-order differential equation, which may be solved 
for a given initial condition of {c n (0)}. An initial value of the expansion coefficient c n (0) is 
determined by 

/oo 
<j) n {x)* V G (x, 0) dx (n = to iV m ), (10) 
-CO 

for a given initial squeezed coherent Gaussian wavepacket ^g(x,0) expressed by [TT1 [T2] 

(ii) 



* G(x ' 0)= (2^)Vi eXP 



[l-ia) 2 .po(x-x )' 

[X — Xq) +1 



4/i v y h 

where x and p are initial position and momentum, respectively, and parameters /i and a 
are related with 

(5x 2 ) = /i, + 5p5x) = a. (12) 

B. Spectral method B 

In an alternative spectral method B, the solution of the time-dependent Schrodinger 
equation may be expressed by 

Nm 

*(x,t) = ^a^(x)e-^\ (13) 

where ^! v {x) and E y are eigenfunction and eigenvalue of the stationary state given by 

HV u (x) = E^ u (x) (i/ = to N m ). (14) 



Note that the expansion coefficient a v in Eq. (13) is time independent and it is determined 
by a given initial Gaussian wavepacket 

/CO 
^g(x,0) dx. (15) 
-co 



Both spectral methods A and B yield the same result. Calculations of time-dependent 



auto-correlation function and various averages expressed by Eq. (30)- (35) in the spectral 
method A are much easier than those in the spectral method B which requires time- 
consuming calculations. By using mostly spectral method A, we have performed model 
calculations to be reported in Sec. II C. The spectral method B is employed for a discussion 
on the tunneling probability in Sec. III. Wavefunctions obtained by the spectral meth- 
ods have been cross-checked by the MATHEMATICA resolver for the partial differential 
equation. 



C. Model calculations 



1. Quartic DW potential 



We consider the quartic DW potential given by jTQl [[3] 

U(x) = (^0 (x 2 - xlf ~d(j- afr) , (16) 

where d expresses a degree of the asymmetry. The potential has locally stable minima at 
x = ±x s and an unstable maximum at x u = d(2x 2 /muj 2 ) with 

U(±x.) = ±^p, (17) 

TT/ \ muj 2 x 2 d 2 x 2 2d^x\ , . 

U(x u ) = -— ^ + § - — ^, 18 

4r/r 3 

AU = u(x.)-U(-x a ) = -^. (19) 
The asymmetry parameter d is assumed to be given by 

7T10J 

-d c <d<d c = — , (20) 
zx s 

for which x u locates at — x s < x u < x s . As for Uo(x), we adopt an HO potential given by 

, x muj 2 x 2 , . 

Uo(x) = — 2— , (21) 

where u) denotes the oscillator frequency. Eigenfunction and eigenvalue for H are given by 



, / n 1 /muj\ 1 / 4 ( mujx 2 \ ( Imuj \ , . 

Mx) = ^(^n) ^(-^rj^-Wx^' (22) 

E 0n =(n+^)hw (n = 0,l,2-,---), (23) 
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FIG. 1: (Color online) (a) The asymmetric DW potential for d = 0.0 (solid curve), d = —0.01 
(dashed curve), d = —0.033 (chain curve) and d = —0.066 (double-chain curve) with x s = 2y/2 
and A = 1.0 in Eq. ( |l6| ). (b) Eigenvalues of E v {y = to 4) as a function of ±d. E v and U{x) are 
symmetric and anti-symmetric, respectively, with respect to a sign of d. 

where H n {x) stands for the Hermite polynomials. We have adopted parameters of m = uj = 1 
and x s = 2^2, which yields d c = 0.1768 and A = U"(x s ) = U^{x s ) = 1.0 for d = 0.0. 
Matrix elements H mn in Eq. ^ are given as follows: We rewrite the potential U(x) as 

U(x) = ^ + ^f + ^f + A lX + A , (24) 

with 

A mUj2 A 7 A mUj2 A 7 2 mU) 2 d 2 

4 = ~2JP' M = -d, A 2 = g", A 1 = dx 2 s , A = ^^. (25) 

After some manipulations with the use of relations given by 

2 V + a), p = i-^=(J -a), (g=±) (26) 



a f (j) n = Vn + 1 (j) n+ i, a (f) n = y/n (f> n -i, (27) 



7 



we obtain the symmetric matrix elements H mn for m < n given by 



3A 4 g 



n + l/2)hcu + —^(2n 2 + 2n + 1) + -^(n + 1/2) + A 



16 



+ 
+ 



1/2 



J m,n— 1 



^-(n - l)^n(n - 1) + - 1 

o 4 



J m,n—2 



+ ^ (|) 7 \/n(n - - 2 ) #m,n-3 

H V n ( n - l )( n - 2 )( n - 3 ) <W-4, 



(28) 



where A' 2 = A 2 — muj 2 . 

Various time-dependent quantities may be expressed in terms of {c n (t)} as follows: After 



some manipulations with the use of the relations Eqs.(26) and (|27|), the auto-correlation 
function is given by 



C(t) 



V(x,t)*V(x,0) dx, 



N„ 



= Y,c n (t)* c n (0). 

n=0 

Expectation values of x(t) and p(t) are expressed by 

ffi? x-^ r 

(P(*)> = W^2][\ / ^+l4 + lWCnW-V^<-lWc„(t) 



(29) 
(30) 



(31) 
(32) 



(^) 2 ) = (|) + + 2) c* n+2 (t)c n (t) + (2n + 1) c* n (t)c n (t) 

n 

+ Vn(n-l)<_ 2 (t)c n (t)], (33) 
Y^W(n + l)(n + 2)c* n+2 (t)c n (t) - (2n + 1) <(t)c n (t) 



2^ 



+ v / n(n-l)c;_ 2 (t)c n (t)], 
(ar(*)p(a;)+p(*)a;(t)> = i fc^h/(n + l)(n + 2) < +2 (t)c n (t) 



^(n-l)<- 2 W)]- 



(34) 



(35) 



The DW potential given by Eq. ( |16[ ) for typical values of d = 0.0 (solid curve), d = —0.01 
(dashed curve) and d = —0.033 (chain curve) is plotted in Fig. [jja). Eigenvalues E v for 



v = — 4 are plotted as a function of ±d in Fig. [jjb). Table 1 shows £7(±x s ), U(x u ), AU 
(= U(x s ) — U(—x s )), S (= Ei — Eq) and 8 f (= i?2 — E\) as a function of an asymmetry d. 
With increasing both \AU\ and 5 are increased. For \d\ > 0.02, 5' becomes smaller than 
8. For d > and d < 0, At/ become AU > and AU < 0, respectively. By using the 
method mentioned in Sec. II, we have studied dynamics of Gaussian wavepackets in DW 
systems. 



d 


U(-x s ) 


U(x u ) 


U(x s ) 


AU 


5 


5' 


0.0 


0.0 


1.0000 


0.0 


0.0 


0.023923 


0.61849 


±0.010 


T0.150849 


1.0064 


±0.150849 


±0.301699 


0.264823 


0.47753 


±0.020 


T0.301699 


1.02555 


±0.301699 


±0.603398 


0.525453 


0.29749 


±0.033 


T0.497803 


1.06929 


±0.497803 


±0.995606 


0.829368 


0.10581 


±0.040 


T0.603398 


1.10153 


±0.603398 


±1.2068 


0.902893 


0.17417 


±0.050 


T0.754247 


1.15787 


±0.754247 


±1.50849 


0.955426 


0.36639 


±0.066 


T0.995606 


1.27231 


±0.995606 


±1.99121 


1.01905 


0.66749 



Table 1 Potential values at locally-stable minima (±x s ) and an unstable maximum po- 
sition (x u ), AU [= U(x s ) — U(—x s )], and energy gaps (5 — E\ — Eq, 5' = E 2 — E\) as a 
function of the asymmetry d for the asymmetric DW potential [Eq. ( fl6| )] (iV m = 30). 



2. Symmetric case 

First we consider the case of the symmetric potential with d — 0.0. A diagonalization of 
the energy matrix with N m = 30 leads to eigenvalues of E v = 0.450203, 0.474126, 1.09262, 
1.39334 and 1.91286 for v = to 4, respectively, which are plotted in Fg. [jjb). The ground 
and first-excited states are quasi-degenerated with a energy gap of 5 = 0.023923. In order 
to examine effects of an initial Gaussian wavepacket on dynamical properties, we consider 



the two Gaussian wavepackets given by Eq. (11) with /i = 0.1 and /i = 0.5 for x = —x S) 



Po = 0.0 and a = 0.0, which are plotted by solid and chain curves, respectively, in Fig. [2} 
The dashed curve will be explained later (Sec. III). 

Figures[3](a) and|3](b) show 3D plots of t)\ 2 calculated by initial Gaussian wavepack- 
ets with ji = 0.1 and /i = 0.5, respectively. As time is developing, initial Gaussian wavepack- 
ets are deformed, and wavepackets at t > cannot be expressed by a single Gaussian [13J. 
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FIG. 2: (Color online) Magnitudes of Gaussian wavepackets (GWs) with n — 0.1 (solid curve) 
and fi = 0.5 (chain curve), dashed curve expressing the two-level wavepacket (TW): \I/(x,0) = 

[^oW-^iW]/^. 

Figures [l](a) and[I](b) show the auto-correlation |C(i)| 2 for initial wavepackets with /x = 
0.1 and \i = 0.5, respectively. Both auto-correlation functions oscillate with a period of 
about 260, which is consistent with the period given by T = 2tt/5 = 262. 

Figures [5]( a) and[5](b) show expectation values of (x) for initial wavepackets with \i — 0.1 
and ji = 0.5, respectively. It is clearly shown that a particle tunnels between the left and 
right wells with a period of about 260. 

Figures |6](a) and[6](b) show the uncertainty product, (5x 2 )(5p 12 ) ) for initial wavepackets 
with n = 0.1 and \i = 0.5, respectively, where 5x = x — (x) and Sp = p — (p) . They start 
from the minimum uncertainty of /i 2 /4 at t = and oscilate with fairly large magnitudes 
and with a period of about 130, a half of the period of (x) in Fig. [5| Its magnitude for 
ji = 0.1 is larger than that for ji = 0.5 by a factor of about four. 



3. Asymmetric case 

Next we consider the asymmetric case with d ^ 0. For d = — 0.01, the potential minimum 
in the right well is lower than that in the left well by AU = —0.301699 (Table 1). We 
obtain eigenvalues of E v = 0.328786, 0.59361, 1.07114, 1.40643 and 1.91312 for v = - 4, 
respectively, which are plotted in Fig. [11(b). Quasi-degeneracy between Eq and E\ for d = 0.0 
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FIG. 3: (Color online) 3D plots of |\I/(x, t) | 2 as functions of x and t for initial Gaussian wavepackets 
with (a) // = 0.1 and (b) \i — 0.5 in the symmetric DW system (xq — — 2y/2,po — 0.0). 

is removed by an introduced asymmetry, while E 2) E% and E± are almost independent of d. 

Figures [7]( a) and [7](b) show 3D plots of |^(x,t)| 2 for d = —0.01 calculated by initial 
Gaussian wavepackets with \i — 0.1 and \i = 0.5, respectively, for x = p G = 0.0 

and a = 0.0. A comparison between Fig. [7](a) [Fig. [7|(b)] and Fig. [3](a) [Fig. [3](b)] shows 
that |^(x,t)| 2 for d = —0.01 is stayed in the left well and tunneling of a particle is almost 
vanishing. This is more clearly seen in Figs. |8](a) and|8](b) which show time dependences of 
(x) for ji = 0.1 and \± = 0.5, respectively.. 

We furthermore increase the asymmetry to d = —0.033, for which eigenvalues are E v = 
0.0193182, 0.848686, 0.954496, 1.46802 and 1.92407, respectively [Fig. [TJb)]. The energy 
gap of 5 (= E x - E ) = 0.829368 is larger than 8' (= E 2 - E x ) = 0.10581, and the 
difference between two potential minima becomes AU = —0.995606 (Table 1). We note 
that AU ~ lj = 1.0, for which a resonance of tunneling is expected. Indeed, expectation 
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FIG. 4: (Color online) The time dependence of the auto-correlation function |C(t)| 2 in the 
symmetric DW system calculated for initial Gaussian wavepackets with (a) n — 0.1 and (b) n — 0.5 
(x = -2V2,po = 0.0). 
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FIG. 5: (Color online) The time dependence of (x) in the symmetric DW system calculated for 
initial Gaussian wavepackets with (a) \x = 0.1 and (b) /j, = 0.5 (xq = — 2\/2, Po = 0.0). 
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FIG. 6: (Color online) The time dependence of the uncertainty product: (8x 2 )(8p 2 ) in the sym- 
metric DW system calculated for initial Gaussian wavepackets with (a) fi = 0.1 and (b) \i — 0.5 
(x = -2 v / 2,Po = 0.0). 



values of (x) for /i = 0.1 (Fig. [8](c)) and ji = 0.5 (Fig. [8](d)) show tunneling with a period 
of about 60. This figure agrees with 2n/8 f = 59.382, which implies that contributions from 
the first- and second-excited states play important roles in the case of d = —0.033. 

Figures [9](a)-[9](d) show time dependences of the uncertainty product of (Sx 2 )(Sp 2 ) for 
d = —0.01 and —0.033, which should be compared to those for d = 0.0 shown in Fig. [6| 
The uncertainty product for d = —0.01 is smaller than that for d = 0.0. It is, however, 
again increased for d = —0.033. We note that (8x 2 )(Sp 2 ) in the resonant tunneling state 
with d = 0.0 or d = —0.033 is larger than that in the non-resonant tunneling state with 
d = —0.01. This is mainly due to the fact that (8x 2 ) in the former state is larger than that 
in the latter. Magnitudes of uncertainty product for a Gaussian wavepacket with fi — 0.1 
are larger than that with \i — 0.5. 

The tunneling probability of P r (t) for finding a particle in the right well is defined by 



p r (t) 



V{x,t)*V(x,t) dt, 



and its maximum by 



P™ ax = max PJt). 



(36) 



(37) 
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FIG. 7: (Color online) 3D plots of |\I/(x,t)| 2 as functions of x and t calculated by Gaussian 
wavepackets with (a) \± — 0.1 and (b) \± — 0.5 in the asymmetric DW system with d = —0.01 
(x = -2V2,p = 0.0). 



Figures flQ|(a), 10(b) and[K)](c) show P r (t) for d = 0.0, —0.01 and —0.033, respectively, which 
are calculated by the initial Gaussian wavepacket with \± — 0.1. For d — 0.0, P r (t) oscillates 
with a period of about 260, as shown in Figs. [3]and[5j For d = —0.01, P r (t) almost stay 
at about 0.2 where it significantly fluctuates. For d = —0.033, P r (t) again oscillates with a 
period of about 60, as shown in Figs. [8](c) and[8]( 



The maximum value of P™ ax is plotted as a function of AU in Fig. 11 where solid and 
dashed curves show the results calculated by initial Gaussian wavepackets with fi — 0.1 and 
\i = 0.5, respectively. The maximum value of P™ ax ~ 1.0 for symmetric DW case [d = 0.0) 
is considerably reduced by an introduced small asymmetry. For a negative AU = —0.995 
(d = —0.033), P r max shows an enhanced value by a resonance effect, while there is no 
resonance for a positive AU = 0.995 (d = 0.033). Similarly, the resonant tunneling is realized 
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FIG. 8: (Color online) The time dependence of (x) in the asymmetric DW system calculated by 
initial Gaussian wavepackets with (a) ji = 0.1 and (b) fi = 0.5 for d — —0.01, and with (c) fi = 0.1 
and (d) n 0.5 for d -0.033. 

for a negative AU = -1.991 (d = -0.066) but for a positive AU = 1.991 (d = 0.066). The 
reduction of P™ ax by the asymmetry for an initial wavepacket with \i = 0.5 is more significant 
than that with \i = 0.1. The AU dependence of P™ ax is not symmetric with respect to a 
sign of AU , which is in contrast with the result of Ref. [5]. 

III. DISCUSSION 

A. The potential-asymmetry dependence of P™ ax 



We will discuss the d (or AU) dependence of the tunneling probability, by using the spec- 



tral method B mentioned in Sec. II B. From Eqs. (13) and (36), the tunneling probability 
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FIG. 9: (Color online) The time dependence of (5x 2 ) (5p 2 ) in the asymmetric DW system calculated 
by initial Gaussian wavepackets with (a) fi = 0.1 and (b) fi = 0.5 for d — —0.01, and with (c) 
fjb = 0.1 and (d) fi = 0.5 for d = -0.033. 

is expressed by 

p r(t) = E E a >A D ^ t/H > (38) 



u=0 A=0 



with 



POO 

D vX = / * u (x)** x (x) dx, (39) 

JO 

where AE„\ = E v — E\. When main contributions arise from the two terms of v — n and 



v = m in Eq. (13), we may adopt the two-level approximation given by 

t) ~ a n * n (x) e - iB »*/ ft + a m * m (a;) e"^' 1 (a* + = 1), (40) 
leading to the tunneling probability 

P r (t) ~ Kl 2 !)™ + |a m | 2 £> mm + 2 ft[a* n a m D nm e ^ E ^l\ (41) 
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FIG. 10: (Color online) The time dependence of P r (t) in the asymmetric DW system for (a) 
d = 0.0, (b) d = —0.01 and (c) d = —0.033 calculated by the initial Gaussian wavepacket with 
fi = 0.1. 



Equations (15), (38) and (39) signify that P r (t) depends on the initial Gaussian wavepacket 



^g( x i 0) and the asymmetry d through the rf-dependent a U) E v and ^ u (x). 

Circles in Figs. [l2|(a) and [l2|(b) show magnitudes of calculated expansion coefficients 
\a„\ 2 of symmetric DW systems (d = 0.0) for initial Gaussian wavepackets with \± — 0.1 



and ji = 0.5, respectively. \a u \ 2 for \i = 0.5 in Fig. 12(b) has main contributions from 
v = and v — 1, which is similar to the conventional two-level wavepacket: 0) = 

ao^o(x) + ai^i(x) with | CLq\ 2 — \cti | 2 = 0.5 plotted by dashed curve in Fig. [2| In contrast, 



\a u \ for /x = 0.1 in Fig. |12|(a) has extra contributions from v = 5 — 7 besides those from 
v = and v = 1 although magnitudes of the former are smaller than those of the latter. 



From Eq. (41) the transition probability for d = 0.0 (AU = 0.0) is given by 



P r (t) ~ al D m + a\ D n + 2a^a x D 01 cos(AE 01 t/h) for AU = 0.0. 



(42) 



Although Eq. (42) leads to a sinusoidal oscillation, P r {t) in Fig. [I0](a) includes fine struc- 
tures, which arise from high-frequency contributions neglected in the two-level approxima- 



tion in Eq. (40). The essential feature of P r (t) for d = 0.0 in Fig. flOfa) may be explained 



by Eq. (42) with \AE 01 \ = S = 0.023923 which yields T = 2n/5 = 262. 
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FIG. 11: (Color online) P™ ax as a function of AU calculated by initial Gaussian wavepackets 
with n — 0.1 (solid curve) and fi = 0.5 (dashed curve). 




k 1 1 1 1 1 r 




FIG. 12: (Color online) Magnitudes of expansion coefficients \a y \ 2 [Eq. (13)] against v with 
d = 0.033 (triangles), d = (circles), d = —0.01 (inverted triangles) and d = —0.033 (squares) for 
initial Gaussian wavepackets with (a) n — 0.1 and (b) \± — 0.5. 

When an asymmetry of d = —0.01 (AU = —0.3016) is introduced, the dominant contri- 
bution comes from v — 1 both for \i = 0.1 and /x = 0.5 as shown by inverted triangles in 



Figs. [I2j(a) and 12(b). This suggests that the wavefunction for d = —0.01 may be given by 
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FIG. 13: (Color online) Gaussian wavepacket of ^g(x,0) (G) with fi = 0.5 and eigenfunctions of 
with v = 0, 1 and 2 for (a) rf = -0.033, (b) rf = 0.0 and (c) d = 0.033. 

the one-level state which yields the time-independent tunneling probability given by 

P r (t) ~ a\ D n for AU = -0.3016. (43) 



Our calculation of P r (t) for AU = —0.3016 (d = —0.01) in Fig. 10 (b) shows wiggles, which 
arise from high-energy contributions not taken into account in the one-level approximation. 
When an asymmetry is increased to d = —0.033 (AU = —0.9956), main contributions to 



\a„\ 2 come from v — 1 and v = 2, as shown by squares in Figs. [l2|(a) and 12(b). From Eq. 



(41), we obtain the tunneling probability 

P r (t) ~ a\ D n + a\ D 22 + 2a x a 2 D 12 cos(AE 12 t/h) for AU = -0.9956. (44) 
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Indeed, P r (t) for d = —0.033 in Fig. |lO|(c) oscillates with a period of about 60 which is 
consistent with T = 2tt/5' = 59.382 for \AE 12 \ = 5 f = 0.10581 (Table 1). 

For a negative A £7, the v — contribution to |a^| 2 is completely suppressed as shown 
in Fig. [12| It is, however, not the case for a positive AU where the v — contribution is 



predominant as shown by triangles for d = 0.033 (AU = 0.996) in Fig. 12 The wavefucntion 
is approximately expressed by the single v — state which yields the time-independent 
tunneling probability 

P r (t) ~ al D 00 for AU = 0.9956. (45) 
The result for a positive AU = 0.9956 is in contrast to that for a negative AU = —0.9956 



given by Eq. (44). 

In the following, we will elucidate the difference between the v dependence of \a u \ 2 for 



AU = —0.9956 and AU = 0.9956, which may be understood from Eq. (15) expressed 
in terms of the eigenfunction of ^ y (x) and the initial Gaussian wavepacket of ^g(x,0). 
Eigenfunctions ^ u (x) (v — to 2) for the asymmetry of d — —0.033, d = 0.0 and d = 0.033 
are plotted in Figs. [l3](a) 7 [l3](b) and[l3|(c), respectively, where \&g(x, 0) with /i = 0.5 is also 



shown. Figure 13(b) shows that the ground-state eigenfunction ^o(^) for d = 0.0 has the 
equal magnitude at x — ±x s . In contrast, ^o(^) fo r d = —0.033 at x = x s has a larger 
magnitude than that at x — x g as shown in Fig. [l3|(a). On the other hand, Fig. [l3|(c) 
shows that the situation is reversed for d = 0.033: |^o( — %s)\ 2 > |^o(^s)| 2 - We note in Fig. 
131(a) that magnitudes of a\ and a 2 for d = —0.033 (AU = —0.9956) may be appreciable 
because ^i(x) and ^2(x) overlap with 0). In contrast, Fig. [l3|(c) shows that a\ and 

a 2 for d = 0.033 (A = 0.9956) become very small because ^i(x) and ^2 0*0 have nodes near 
the center of ^g(x,0) while a is appreciable because ^o(^) and ^g(x,0) are overlapped. 
It is necessary to note that high-energy contributions to |a^| 2 at v = 5 — 7 for ji = 0.1 are 



almost independent of the asymmetry in Fig. [121(a) while there are no such high-energy 



max 
r 



contributions for /i = 0.5 in Fig. 12(b). This is the reason why the AU dependence of P 7 
for n = 0.1 is smaller than that for /x = 0.5 as shown in Fig. [TTj 

When the asymmetry is much increased up to d = ±0.066 (A = ±1.991), the energy gap 
between the second- and third-excited states: E 3 — E 2 = 0.230984 becomes smaller than 5 
and 5' with E 2 (= 1.25498) < U(x u ) (= 1.27231) < E 3 (= 1.48597) [see Fig. [I] (b)]. Then 
for a negative A = —1.991 (d = —0.066) dominant contributions to |a^| 2 arise from two 
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levels of v — 2 and v = 3, while for a positive A = 1.991 (d = 0.066) a contribution from 
a single level of v = is predominant (relevant results not shown). This is similar to the 
case of AU = —0.9956 and AU = 0.9956 mentioned above. We may similarly elucidate the 
difference between P™ ax of AU = -1.991 and AU = 1.991 in Fig. [TT, 



Cordes and Das (CD) [8] discussed the tunneling probability in asymmetric DW systems, 
proposing the generalized two-level wavefunction given by 



■iEmt/h 



(a 2 n + a 2 m = 1). 



(46) 



Here eigenfunctions ^^ D (x) and ^^(x) of the DW system are assumed to be expressed 
by superposition of eigenfunctions for two harmonic potentials in left and right wells which 



are separated by a high central barrier. By using Eq. (46), CD showed that the tunneling 
probability is given by j8] 



P r CD (t) = 2a 2 n a 2 m [l - cos(A£ t/h)\ (AE = \E m - E m \). 



(47) 



The tunneling probability given by Eq. (47) is consistent with our results for d = 0.0 and 



d = —0.033 given by Eqs. (42) and (44), respectively. However Eq. (47) is not valid for 
cases of d = — 0.01 and d = 0.033 for which it yields P^ D (t) = in contrast to Eqs. (43) 



and (45). Actually, P^ D (t) in Eq. (47) cannot be applied to the one-level state with either 
a n = or a m = 0, while P r {t) given by Eq. (41) is applicable. 



B. Varying initial Gaussian wavepacket 



In our study having been reported in Sees. II and III, we have adopted the initial squeezed 
Gaussian wavepacket given by Eq. (11) with x = —2^/2, Po = 0.0, a = 0.0 and /x = 0.1 (or 



\i = 0.5). We may, however, employ an any arbitrary initial wavepacket with appropriate 
parameters of x , p , /x and a, while conventional theories rely on the two-level wavepacket. 
For example, we here employ an initial Gaussian wavepacket with x = 0.0, po = 0.5, ji — 0.1 
and a = 0.0. In the classical mechanics, a particle with this initial condition rolls down from 
the origin near a top of the potential with an initial velocity of v = po/m = 0.5, and it 
continues an oscillation at —4 < x < 4. Our calculation, however, shows that motion of the 
particle in quantum mechanics is different: (x) for d = 0.0 and d = —0.033 stay near the 



origin (\(x)\ < 2) as shown in Figs. [L4[a) and 14(b), respectively. We note that results for 
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FIG. 14: (Color online) The time dependence of (x) for (a) d = 0.0 and (b) d = —0.033 calculated 
by an initial Gaussian wavepacket with xo = 0.0, po = 0.5, /i = 0.1 and a = 0.0. 

the initial condition of x = 0.0 and p = 0.5 in Figs. [l4](a) and [l4|(b) are quite different 
from relevant results for x = —2^/2 and p = 0.0 shown in Figs. [H](a) and[8](c). Thus the 
time dependence of (x) and other quantities depend on the assumed initial condition. 



IV. CONCLUSION 



Dynamics of Gaussian wavepackets and quantum tunneling in asymmetric DW systems 
have been studied with the use of the numerical method which has advantages that (a) it is 
simple and physically transparent, (b) it is applicable to realistic DW potentials, and (c) it 
may adopt an arbitrary, appropriate initial state. Our calculations have shown the following: 

(1) The maximum tunneling probability P™ ax is considerably reduced by a small amount of 
the asymmetry AU in the DW potential, 

(2) A resonant tunneling with \AU\ ~ k Huj (k = 0, 1, • • • ) is possible for a motion from 
upper minimum (AU < 0) but for a motion from lower minimum (AU > 0), 

(3) P™ ax for an initial Gaussian wavepacket with narrower width (\i) is less sensitive to the 
asymmetry, and 

(4) The uncertainty product in the resonant tunneling state is larger than that in the non- 
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resonant tunneling state. 

The item (1) for AU < is consistent with results of previous studies [3j[U|8]. The item (2) 
is in contrast with Ref. [5j which claimed the symmetric behavior for a motion starting from 
the upper and lower minima. The item (3) is due to the fact that the Gaussian wavepacket 
with a small ji (= 0.1) includes high-energy contributions to a v whose magnitudes are 



nearly independent of the asymmetry (Fig. 12). The item (4) signifies that tunneling 



and uncertainty, both of which are typical quantum phenomena, are mutually related. In 
order to examine a validity of items (l)-(4), it would be interesting to observe |^(x,t)| 2 
in asymmetric DW systems, which seems difficult but possible with the recent advance of 
experimental methods. The present study has been made without considering dissipative 
effects which are expected to play important roles in stationary and dynamical properties of 
real DW systems. An inclusion of dissipation arising from environments is left as our future 
subject. 
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